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In previous work ( physics/ 0004 026 ) was shown, with the help of numer- 



ical calculations, that the effective Brown temperature for charged particle is 
lower than that for particle without charge. Here we derive this result with- 
out numerical calculations, integrating the delay equation analytically, as for 
zero, so for nonzero viscosity. 
03.50.De 

1. 

To describe motion of charged Brown particle in so called "extended 
quasi- stationary approximation" [1] in [2] was used the Sommerfeld model [3] 
of charged rigid sphere. The equation of straightline motion of such Brown 
particle in dimensionless form reads [2]: 

y(x) = f(x) + 7 • [y(x - 5) - y(x)] (1) 

here 

y(x) - is dimensionless velocity of the particle; 
x - is dimensionless "time"; 
f(x) - is some external (stochastic) force; 
5 - is "time" delay; 

7 - is coefficient: 7 • S is proportional to the ratio of particle's electromag- 
netic mass to the mechanical mass: 7 ■ 5 = (2/3)(Q 2 /a)/(mc 2 ) ( 2a - is the 
size of Sommerfeld particle of charge Q and mass m); 

the viscosity T of the surrounding medium is zero. 

In [2] was shown, with the help of numerical calculations, that the effective 
Brown temperature for charged particle is lower than that for particle without 
charge. Here we derive this result without numerical calculations, integrating 
the delay equation (1) analytically. 
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With zero initial conditions: 
y = y = 0foix<0 

dividing the x-axis into 5 - intervals (i — 1)5 < x < i5, i — 1, 
and integrating eq. (1) step by step with boundary conditions yi(x = 
iS) = Vi+i(x = iS), we finally get the recurrence formula: 
for (N -l)5<x<N5 

y(x) = y N (x) = 

rx rx—5 

\ dzf(z) exp 7(2; — x) + 7 / dz yjv-i(^) ex P l( z + 8 — x ) 

JO J(N-2)S 

dz Vi( z ) ex P 7(2 + S - x ) ( 2 ) 

i=l • / (i-l)'5 



with 



(x) = / dzf(z) exp 7(2; — rr), < x < 5 
Jo 



Let's consider one interesting case: 

f(x) for intervals (i — 1)5 < x < i5 is constant and is equal to /j. 
Then the eq.(2) for x = iV<5 = xn yields 

1 JV AT 

y* N = y N (x = x N ) = - ]T fk [1 - C(iV - fc;p)] = ]T /* (3) 
7 fc=i fc=i 

where the function C (n; p) is defined as 

n 

C(n;p) = exp (—p(n + 1)) ^ (pexpp) m (n+1 — m) m /(m\); (4) 

m=0 

here p = jS. 

Function C(n;p) is positive and for sufficiently large n (for ex., n > 20 
for p — 1.0 ) is equival to yrj^ : 

C(^;p)n»i = -4 ► ~ n -j v (5) 

1+p (I+J9J7 

Thus if fi = / =c o n st Vi, then from (3,5) we get for iV 3> 1 



in accordance with the exact solution of (1) for / = / =const: 

y(x) = — x 

Also for N ^> 1 one can rewrite (3) in the form 

p N 5 N 

This result resembles the classical Brown result: from eq.(4) with 7 = one 
immediately gets 

y(x)= [ X f(z)dz, (7) 
Jo 

dividing x-interval of integration in (7) into 5- intervals with f(x) = for 
(k — 1)5 < x < k5, one can take the integral in (7) in the following manner: 

N 

y(x = x N ) = 5J2 fk (8) 

k=l 

This result differs from (6) only in the multiplier jj^j- 

Thus one can say that the effect of delay (effect of retardation) for eq.(l) 
reduces to the effect of mass renormalization: m — > mj (1+p), or consequently 
to the effect of reduction of the external force: 

/-> (9) 

This result also says that the reduction of the external force is model-independent 
one, and instead of 7$ one can write the classical ratio of self-electromagnetic 
mass to the mechanical mass m in its general form: 

7 ^J_/ dr ^>(f) (io) 
mc z J |r — r'| 

here p - is distribution of charge of a particle. 

If fk, k — 1, is the range of stochastic numbers with average value f a : 
< fk >= f a (here brackets <> denote time average with the same definition 
as in the classical theory of Brownian motion), then eq.(3) yields 

N 

<V*N>= faJ2 D k- faX N /(l+p) (11) 
fc=l 
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Consequently the dispersion D is 

N N 

D = (y* N - <y* N >) 2 = J2Y: D k D m < (f k - f a )(f m - f a ) > 

k=l m=l 
N N 

= Y,Y, D kD m R(k-m) (12) 

k=l m=l 

here R{k — m) - is correlation function of stochastic force /. If R is compact: 

R(k -m)= R 5 mk /5 (13) 

then the dispersion (12) is 

N 

D = Ro/S^Dk) 2 ~ R x n /(1+ P ) 2 (14) 
fc=i 

This result should be compared with classical one. 

The theory of Brownian motion without viscosity tells ( eq. (1) with 
7 = 0) that the dispersion D B is 

rx rx 

D B = / dzi / dz 2 ■ R(zi - z 2 ) (15) 

JO JO 

here R{z\ — z 2 ) =< (f(zi) — f a )(f(z2) — fa) > - is the correlation function. 
If 

R(zi - z 2 ) = Rq5(z 1 - z 2 ) 

then 

D B = R x (16) 

Consequently we see that (eqs. (16) and (14) ) the dispersion of the Som- 
merfeld charged particle is lower than that of the classical Brown particle 
without electric charge: D = D B (l +p)~ 2 - Thus one can say that the ef- 
fective temperature of Sommerfeld particle is lower than that of the Brown 
one. This result is model independent one (see the remark made above - eq. 
(10)). 

So we confirm the result of the work [2] . 
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2. 

If the viscosity T is not zero, the main equation reads: 

y(x) + T • y(x) = f{x) + 7 • W - 5) - y(x)} (17) 
For f — fo — const eq.(17) has the exact solution 

y(x) = ^(l-e X p(-ax)) (18) 

and a is determined by the eq. 

r + 7 -a = 7exp(a5) (19) 

Iterative solution y(xjv) = y* N of eq.(17), if f(x) = f = const for intervals 
(i — 1)5 < x < i5, can be put in the form: 

N 

y*N = H f^k (20) 
k=l 

here D k - some discrete function which can be found from recurrence for- 
mula, analogous to (2). But it is convenient to find D k from the following 
considerations, using exact results (18,19). Solution (20) must tend to the 
exact solution (18) (in the case f — fo — const V i) if the a;- axis is divided 
into infinitesimally small (^-intervals: 5^0 and N — > oo in such a way that 
xn = 8 ■ N = const. Thus one can rewrite 



N N 

^2 fk^k — fo X] 
k=l k=l 



V*n = J2 fkD k = f D k = ^(1 - exp (-ax N )) 



so 



N y 

J2D k = -(l-exp(-a5N)) (21) 

k=l L 

If 5 — > we can replace the sum in lhs of (21) by the integral: 

N r N 1 

J2D k ~J D k dk = -(l-exp(-a5N)) (22) 

k=l 

Differentiation of (22) with respect to iV provides us with this expression for 
D N : 
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Av ~ exp (-aSN) (23) 
Substitution of (23) back into (21) gives 

aS a5 l — exp(—a5N) 
> — exp(-a5k) = (24) 

Consequently the required result (rhs of (21) ) is reproduced if we expand 
the denominator in (24) in the following way: 

exp (aS) - 1 « aS (24) 

Using this representation of Dj,, one can find the dispersion D. For cor- 
relation function (13) we have 

k=i 

here we expanded the expression exp (2a5) — 1 in the same manner as in (24): 
exp (2a5) — Is* 2a5. 

Solving the eq.(19) in approximation (24), we find 

(26) 



(l + 7 5) 

So with (26) and (25) the dispersion D takes the form 

D = 2 ( m + W 1 - e * p( - 2aSN)) (27) 

Dispersion (27) for 7 = is exactly the same as Brownian dispersion D B : 

D B = ^(l-exp(-2Tx N )) 
If (ixn ^ 1) solution (27) yields 

RqXn 



i.e. the solution we have got earlier (14). 
If ax N > 1, (27) yields 



~2r(i + 7 5) (i + 7<5) 

Thus dispersion D differs from the Brownian one. Consequently the effec- 
tive temperature of charged particle, undergoing Brownian motion, is lower 
then that of particle without charge. Now we have proved this result in 
general case of nonzero viscosity. Of course, our general conclusion is model- 
independent one - see the above remark (10). 
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